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ABSTRACT 

The installation of the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) will revolution- 
ize the study of high-redshift galaxy populations. Initial observations of the HST Ultra Deep Field (UDF) have 
yielded multiple z > 7 dropout candidates. Supplemented by the Great Observatories Origins Deep Survey 
(GOODS) Early Release Science (ERS) and further UDF pointings, these data will provide crucial information 
about the most distant known galaxies. However, achieving tight constraints on the z ^ 7 galaxy luminosity 
function (LF) will require even more ambitious photometric surveys. Using a Fisher matrix approach to fully 
account for Poisson and cosmic sample variance, as well as covariances in the data, we estimate the uncer- 
tainties on LF parameters achieved by surveys of a given area and depth. Applying this method to WFC3 
z^l dropout galaxy samples, we forecast the LF parameter uncertainties for a variety of model surveys. We 
demonstrate that performing a wide area 1 deg^) survey to Hab ^ 27 depth or increasing the UDF depth 
to Hab ^ 30 provides excellent constraints on the high-z LF when combined with the existing UDF GO and 
GOODS ERS data. We also show that the shape of the matter power spectrum may limit the possible gain 
of splitting wide area (> 0.5 deg^) high-redshift surveys into multiple fields to probe statistically independent 
regions; the increased root-mean-squared density fluctuations in smaller volumes mostly offset the improved 
variance gained from independent samples. 

Subject headings: surveys-methods: statistical-galaxies:abundances 



1. INTRODUCTION 

Recent studies using Hubble Space Telescope (HST) Wide 
Field Camera 3 (WFC3) observations have discovered tens 
of candidate galaxies at redshifts z > 7 ( Bouwens et al] 
20091 I2010 allbf 'Oe schetal.1 l2010b lal: ^Buriker et_aL . 20091 
McLure et a l. 2010: lYanetal.1 1200 9: Wilkins etMl l2010t 
Labbe et alj i2009i, i2010t iFinkelstein et al.i |2009 ) . The new 
WFC3 observations have broadened our knowledge of the 
highest redshift galaxies yet found, complementing z > 
7 galaxy searches with the Near Infrared Camera and 
Multi - Obiec t Spectromete r ( Kneib et al. 2004; Bouwens et al. 
IpP, 12005"; 'Eaami et al. 2005; Henrv et al. 2007, 2008, 
2009; Richard et al. 2008; Bradley et al. 2008; Bouwens et al. 
2008t iZheng et all 120091: lOesch et alj l2009; Gonzalez et al. 
20101) . ground-based dro pout sele ctions (Richard et al. 2006 ; 
Stanwavetal. 2008; Ouc hi etalJ l2009a: H ickey et al.1 bo'iol: 
Castellano et al. 2010), and narrow band Lyman-a (Lya) 
emission s urveys (Parkesetal. 1994; Kodairaetal. 

2004; Wiflis & Courbin 2005; WilHs et al 
2005; Stark & Ellis 2006; lye et al. 
. . , 2007a; Cuby et al 

lOuchi et al 



Santos et al. 
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Kashikawa et all 12006 : Stark et al. 
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Hibon et al. 200^: 



.2008t lOuchi et aL ,2009b, 

Sobral et al.l l2009h . The existence of star-forming galaxies at 
z > 7 has been well-established by these studies, and the im- 
portance of these high-redshift galaxies for reionization and 
subsequent galaxy formation at lower redshifts will likely 
motivate the dedication of large telescope allocations to de- 
tailing their abundance. The purpose of this paper is to de- 
velop a method to rapidly compare possible photometric sur- 
vey strategies for detecting large numbers of z ^ 7 galaxies, 
and forecast constraints on the z ^ 7 galaxy luminosity func- 
tion (LF) achieved by different HST survey designs. 
A determination of the constraining power a survey can 
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obtain requires an estimate of the uncertainty in the abun- 
dance of galaxies as a function of luminosity. In addition 
to the Poisson uncertainty inherent in galaxy counts, cosmic 
sample variance induced by density fluctuations and galaxy 
clustering must b e accounted for (see, e.g., Newman & Davii 
I2OOI iSomervilie et allbOOl fStark et al.ll2007bl) . A partic- 
ularly powerful approach for estimating these uncertainties 
and determining the resulting potential constraints on the lu- 
minosity function of galaxies in the UDF was presented by 
iTrenti & Stiavelli (2008). These authors used cosmological 
simulations to determine the abundance and spatial distribu- 
tion of dark matter halos and then applied a model for the 
halo mass-to-light ratio to determine the abundance of galax- 
ies of a given luminosity. Poisson and sample variance un- 
certainties were estimated by drawing pencil beam realiza- 
tions of the survey from the cosmological volume (see also 
Hitzbichler & White 2007). A maximum-likelihood approach 
was then used to study constraints on the luminosity function 
for various dropout selections. 

The maximum-likelihood estimation of LF parameters 
based on mock catalogues can account for detailed selection 
effects and spatial correlations in addition to the Poisson and 
sample variance uncertainties. Such simulations have clear 
advantages for estimating the completeness of magnitude- 
limited surveys or understanding systematic effects intro- 
duced by dropout color selections. However, the need for 
halo catalogues from cosmological simulations introduces 
two limitations. First, the volume of the simulation should 
probe many independent realizations of the modeled survey. 
However, future high-redshift galaxy surveys with WFC3 
or other instruments may take the form of deeper versions 
of surveys like t he Spitzer Extended Deep Survey (SEDS; 
iFazio et al.l2008l) Exploration Science program or the Cosmic 
Evolution Survey (COSMOS; Scovifle et al. 2007b a). These 
surveys are each ~ 1 deg^ or larger, and large volume (L > 
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200/i 'Mpc) cosmological simulations are required to probe 
multiple independent samples of the surveys' high-redshift 
galaxy populations. For instance, a ^ 1 deg^ survey at 6.5 < 
z < 7.5 has a comoving vol ume of y 3 X 10' ' /r^Mp cl The 
largest simulation used by iTrenti & Stiavellil ( 120081) would 
provide less than two independent samp les of such a volume , 
and even the Millennium Simulation (iSpringel et alJ l2005h 
with a comoving box size L = 500/;"' Mpc would only pro- 
vide ^ 40 independent samples of such a wide high-redshift 
survey. Second, the method requires access to and manipu- 
lation of cosmological simulation results. This requirement 
may pose an unwanted computational overhead for those in- 
terested in rapid estimates and comparisons of potential con- 
straints from a wide range of survey designs. 

A simpler methodology for estimating survey constraints 
on the abundance of high-redshift galaxies that does not di- 
rectly require halo catalogues from cosmological simulations 
is therefore desirable for performing rapid comparisons of 
survey designs. Hence, we seek an approximate method 
for forecasting LF parameter constraints that relies on de- 
scriptions of galaxy and dark matter halo abundance and 
clustering, Poisson and cosmic sample variance, and pa- 
rameter covariances that are analytical or easily calculable 
through numerical methods. We utilize a simple model for 
describing the clustering of z ^ 7 galaxies based on fidu- 
cial empirical estimates of the hig h-redshift galaxy luminos- 
ity function jOesch et al.l 12010a) and abunda nce matching 
between galaxies and dark mat ter halos (e.g., I Conrov et alj 
120061: IConrov & Wechsleiil2009h . We then adopt a common 
approach to translate galaxy clustering and matter fluctua- 
tions into an estimate of the cos mic sample variance (see 
the various cal culations in, e.g., 'Newman & Davis' '2002'; 
Somerville et al.,,200 4; Stark et al.l ,2007b; Tren ti & S tiaveUi, 
20081) . With this estimate of the sample variance and Poisson 
uncertainty from an assumed fiducial model for the abundance 
of galaxies, we use a Fisher matrix formalism to characterize 
the likelihood function and estimate z^l luminosity function 
parameter constraints. The presented methodology is fast and 
flexible, and can be used with appropriate extensions, to esti- 
mate constraints on galaxy abundance for other survey sample 
selections and redshifts. 

Motivated by the exciting recent HST WFC3 results, we 
focus on modeling z^l dropout survey designs. While we 
choose to study broadband searches for high-redshift galaxies, 
narrow band surveys for high-redshift Lya emission present 
another interesting class of survey designs. The rapid progress 
in detecting increasing numbers of high-redshift Lya emitters 
using narrow band surveys has motivated theoretical efforts 
both to understand and predict the abundance of Lya emit- 
ters. The observable properties of the high-redshift Lya emit- 
ter population are particularly difficult to model owing to the 
uncertain escape fraction, intergalactic medium absorption 
and other radiative transfer effects, as well as uncertai nties in 
our knowledge of the galaxy formation process (e.g., Haiman' 
2002; Santos 2004; Barton et al. 2004; Wvithe & Loeb 
2005; Le Defliou et al. 2006; Hansen & Oh 2006; Dave et al. 
2006; Furianetto et al. 2006; Ta sitsiomi 200 6: McOuinn et al. 
20071 INilsson et al. 2007; M ao et al.l 120071: i Kobavashi et al. 
2007", '2010; Stark et al. 2007b; 'Mesinger & Furianetto 200& 
Fernandez & Komatsu 2008; Tilvi et al. 2009; Daval et at] 
2009L I2010t iSamui et al.l 12009 ). While the astrophysics in- 
volved in these studies are tremendously interesting and pro- 
vide another route to probe high-redshift galaxies, we will 



only examine the statistical constraining power of various 
broadband survey designs and will not attempt to model the 
Lya emitter population. 

This paper is organized as follows. Forecasting constraints 
on z ^ 7 LF function parameters requires a model of the 
sources of error and covariances in the data. In ^ we discuss 
sample variance uncertainties owing cosmic density fluctua- 
tions. In ^ we review the Fisher matrix formalism and show 
how to apply the formalism to forecast LF parameter uncer- 
tainties accounting for cosmic sample and Poisson variances. 
To perform actual forecasts for the z ^ 7 LF, we review exist- 
ing HST WFC3 survey data in ^ and define fiducial model 
surveys in ^ In ^ we combine the expected constraints 
from existing surveys with forecasts of luminosity function 
constraints from model surveys. We discuss our results and 
possible caveats in ^ and summarize and conclude in ^ 

Throughout, we work in the context of a ACDM cos- 
mology consistent with joint constraints from the 5-year 
Wilkinson Microwave Anisotropy Probe, Type la supernovae, 
B aryo n Acoustic Osci ll ation, and Hubble Key Project data 
(Free dman et al.1 120011; iPercival et all 120071; iKowalski et all 
l2008l) . Specifically, we adopt a Hubble parameter h = 0.705, 
matter density 0„i = 0.274, dark energy density JIa = 0.726, 
baryon density ilb = 0.0456, relativistic species density il, = 
4.15 X 10"^, spectral index = 0.96, and root-mean-squared 
density fluctuations i n 8 /!"'Mpc-radius spheres of (j% = 0.812 
(|Komatsu et al. 2009). We report all magnitudes in the AB 
system (Oke & Gunn 1983). 

2. POISSON UNCERTAINTY, COSMIC SAMPLE VARIANCE, AND 
THE ACDM POWER SPECTRUM 

We wish to evaluate the relative merits of various galaxy 
survey designs in terms of their ability to constrain the galaxy 
luminosity function. To perform this evaluation, we must de- 
termine the quality of each design in terms of the number of 
galaxies of a given luminosity the survey will discover (the 
Poisson variance) and the intrinsic scatter expected for the 
survey volume given variations in the cosmological density 
field (the cosmic sample variance). This section of the paper 
formally defines each source of uncertainty and describes how 
these variances are calculated. 

We define cosmic sample variance as the fluctuations in a 
volume-averaged quantity owing to density inhomogeneities 
seeded by the matter power spectrum. We will use the terms 
"cosmic variance" and "sample variance" interchangeably, 
but elsewhere in the literature cosmic variance is taken to 
equal the sample varianc e only in the lirnit of the entire vol- 
ume of the universe (e.g.. lHu & Kravtsovll2003h . 

2.1. Dark Matter Density Variance 

Density fluctuations, or differences between the local mat- 
ter density p,„(x) and the mean matter density p^, can be de- 
scribed in terms of a local matter overdensity 5(x) = [/9m(x)- 
Pm]/Pm- Consider a survey of comoving volume V at redshift 
z. For "unbiased" quantities measured within V that spatially 
cluster like the dark matter, such that the two-point correlation 
function ^(r) is identical to the dark matter correlation func- 
tion ^,„(r) = {S{x)S(x + r)), the cosmic variance is simply the 
dark matter variance 

D\z)<7l^={S^(x))=D\z) I -^P(A:)|W(k,y)|2, (1) 

where W'(k,V) is the Fourier transform of the survey vol- 
ume window function W(x) (whose geometry may introduce 
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a dependence on the direction of the wavenumber k), D{z) 
is the linear growth function, and P(k) is the isotropic lin- 
ear ACDM power spectr um. To c alculate P(k), we use the 
transfer function of Eisen stein & H u (1998) that includes the 
effects of baryons. We ignore possible nonlinear correc- 
tions to the powe r spectrum (e.g.. Peacock & Dodds 1996; 
ISmith et ani2003h . The window function is normalized such 
that / d?xW{x) = 1 . For a spherical volume of comoving ra- 
dius R = 8/!"'Mpc, Equation[T] would provide ctdm = red- 
shift z = 0. The linear growth function 

/■°° (l+z')dz' 
D(z) = D,mz)j '-j^ (2) 

has a normalization constant Dq chosen such that D(z = 0) = 1 . 
The Hubble parameter 

H(z)=Ho[n,(i+z)^+nM+z)' 

+(i-r!,„-r!A-a)(i+z)'+f^A]'/' (3) 

describes the rate of change of the universal scale factor 
// = fl'/fl as a function of the matter density fl^, relativistic 
species density fti-, and dark energy density (taken to be a 
cosmological constant). 

2.2. Dark Matter Halo Abundance and Clustering 

For galaxy surveys, where quantities of interest depend on 
the abundance and clustering of galaxies, the sample variance 
will depend on the bias of dark matter halos hosting the ob- 
served systems. We can define the bias in terms of the corre- 
lation function as = ^/,/^,„ where is the correlation func- 
tion of dark matter halos. Given a halo mass function, the 
bias of halos with mass m can be estimated using the peak- 
background split formalism (e.g., Kaiserl ll984l: iMo & White! 
11996; Sheth & Tor men 1999) or measured directly from the 
simulations via the halo correlation function or halo power 
spectrum. We adopt the latter approach. 

We use the dark matter halo mass function measured by 
iTinker et al.' (W0§) fr om a large suite of cosmologica l simu- 
lations (Kravtsov et al. 2004; Warren et al. 2006; Croc ce et al.l 
2006; Gottlober & Yepesj 120071: lYepes et all 120071) . The 
Tinker et al. (2008) mass function can be written as a function 
of the dark matter halo mass m in terms of the "peak height", 

D(z)a(m) 

where Sc = 1 . 686 is the spherical collapse barrier (see, e.g, 
iGunn & GottI [T97l IB ond & MversI [19961) . D(zMm) is the 
square root of the dark matter variance (Equation[T]i evaluated 
in a spherical volume of comoving radius R = (3m/47r 
Here, p„, is the background matter density. The linear growth 
function D{z) is given by Equation|2] 

With the definit ion of peak height ly in Equation H] the 
ITinker et alj (l2008h halo mass function can be written 

-i- = —f(^):^^ (5) 

am m am 

where the function 

/(z/) = a[l+(/3j.)-2^]i^2^e-^'''/2, (6) 

is ofte n called the "first crossing distribution." ITinker et al.l 
( 120081) find that the parameter values a = 0.245, /3 = 0.757, 
7 = 0.853, = -0.659, and r/ = -0.341 fit well the z = 2.5 



simulated mass function measured for halos defined with a 
spherical overdensity A = 20 relative to the b ackground mat- 
ter density (see also §4 of ITinker et al.l l2010i) . The abun- 
dance of dark matter halos more massive than m is then just 
nh( > m) = (d«h / dm)dOT. 

We choose the Tinker et"an (120081) mass function because it 
is accurate to < 5% for halos in the mass range lO"/z~'M0 < 
m < 10'''/!~'Mq at redshift z = 0, and improves on previ- 
ous approximations by 10-20% (c.f., Sheth & Tormen 1999). 
ITinker et al.l (120081) demonstrate that the halo mass function 
does not have a redshift-independent, universal form, and that 
the normalization of the first crossing distribution evolves at 
the 20-50% level between z = Q and z = 2.5. However, 
the halo mass function has not been calibrated at the red- 
shifts of interest (z ^ 2.5) and following the advice in §4 of 
[tinker et al. (2008) we will use the z = 2.5 first crossing distri- 
bution as the best available approximationO We note that us- 
ing any other previously published m ass function from the lit- 
erature (e.g.. lSheth & Tormenll 19991) will therefore introduce 
an unknown error in the abundance of halos at high redshifts. 
Tinker etal. (2008) estimate this error could be as large as 
^ 20-50% for galaxy-sized halosQ 

For the bias b relati ng halo and dark ma tter clustering, we 
will use the results of ITinker et al.l ( 120101) who measure the 
halo bias as a function of peak heig ht in a m anner consistent 
with the halo mass function of Tin ker et"an (|2008). The bias 
function biiy) is constrained by the halo first crossing distri- 
bution f( i') by requiring that dark matter is not biased against 
itself, i.e., 

j b{v)f{u)Ay=\. (7) 

Under this constraint, ITinker et al.l (120101) find that the fitting 
function 

b{v)=\-A^^+Bv'' + Cv' (8) 

u" + 0° 

with parameters A = 1.0, a = 0.1325, B = 0.183, b= 1.5, 
C = 0.265, and c = 2.4 provides an accurate match to the 
bias of dark matter halos defined with a spherical overden- 
sity of A = 200 rel ative to the b ackground matter density. 
As demonstrated by I Tinker et~al] {2010), Equation |8] repro- 
duces the simulated halo clustering better than the analytical 
formulae of Mo & White ( 1996) or Sheth et al. (2001) calcu- 
lated u sing the peak-background split formalism. ITinker et al.l 
(120101) find that the bias b{v) as a function of peak height v is 
nearly redshift-independent, and we will adopt Equation|8]for 
b{i/) at all redshifts. 

2.3. Galaxy Abundance and Clustering 

Our main premise is to use the Fisher matrix approach to es- 
timate the constraints on a model for the abundance of galax- 
ies that reproduces well the observed source counts. The def- 
inition of the model for galaxy abundance is therefore im- 

^ While we account for the redshift-dependent abundance of dark matter 
halos at the mean redshift of the survey, we ignore the evolution of the dark 
mat ter halo abundan ce over the redshift interval of the survey volume. See, 
e.g. JMuiioz & Loebl (2008) for a quantification of the effect of this evolution 
on the variance of inferred halo number densities. 

^ We find that when using the Sheth & Tormen 1 1999) mass function and 
the corresponding Sheth et al. (2001) bias function, the marginalized errors 
calculated in sections and |7] degrade by ~ 5% relative to the results ob- 
tained with the Tinker et al. 12008) mass function and Tinker et al. (2010|) 
bias function. This difference quantifies how the results of the presented 
method depend on the choice for the halo mass and bias functions. 



4 



portant. Further, the clustering of galaxies directly influences 
the covariances of the data and knowledge of the galaxy spa- 
tial distribution is therefore also quite important. Since we 
are interested in the characteristics of the high-redshift galaxy 
population for which little clustering information is known, 
we will estimate the galaxy clustering bias by associating lu- 
minous galaxies with dark matter halos of similar comoving 
abundance and assigning those galaxies the bias of their asso- 
ciated halos (as provided by Equation [8]). When more de- 
tailed clustering information is available, as is the case at 
lower redshifts, additional constraints on the connection be- 
tween gal axy and halo populations are attainable (see, e.g., 

iLee et al.ll2009.) . 

We will adopt the commonlv-used lSchechterl (119761) model 
for the abundance of galaxies. Specifically, the expected num- 
ber density n, of galaxies in the i-th luminosity bin of width 
AM about magnitude M, can be written 



j.Mi+AM/2 

ni= / $(M)dM, 

Jmi-AM/2 



where the iSchechteii i 1 9761) function 



$(M)=-ln(lO)0, 



2Q§(M.-M) 



a+l 



exp 



-105 



(M.-M) 



(9) 



,(10) 



describes the distribution of galaxy luminosities, is the lu- 
minosity function normalization in comoving Mpc~^ mag"', 
is the characteristic galaxy luminosity in AB magnitudes, 
and a is the faint-end slope. W e will often refer to the param- 
eters of this ISchechterl ( 119761) model in terms of the vector 
p = [logiQ(j)i,,Mi,,a], and it is these parameters for which we 
will forecast constraints. The fiducial values for the parame- 
ters p used in the Fisher matrix calculation will be selected in 

m 

In a manner similar to Equation|9] we can also define the co- 
moving abundance Hi of galaxies more luminous than magni- 
tude M as «/,(< M) = $(M)dM, where the negative lower 
limit follows from the definition of magnitudes. 

With a model for the abundance of galaxies, we will as- 
sociate galaxies with dark matter halos of similar abundance 
to estimate the galaxies' spatial clustering. The abundance 
rii of galaxies in the range M, ± AM/2 can be written as 
n,- = «/,(< Mi + AM/2) - niiK Mi - AM/2). We match the 
abundance of galaxies and halos as 



nh(> m) ~ nL(<M) 



(11) 



at the minimum and maximum luminosity of galaxies in each 
magnitude bin (e.g., M = M,- + AM /2 and M = M; - AM /2), 
which provides a mas s range ot, zfc Aw/2 of halos with a sim- 
ilar abundance (e.g., Conrov et al. 2006^ IConroy & WechsleJ 
12009.) . The comoving number density of these halos is sim- 
ply «h,! = «h(> Mi — Am/2) — nh(> m, + Aot/2), with nh,; = «!■ 
The resulting connection between galaxy luminosity and halo 
mass is simplistic, but more sophisticated stellar mass-halo 
mass relations could be incorporated into our approach when 
warranted by th e constraining power of the available data 
(see, e.g., Behr oozi et al.ll201d) . We adopt AM = 0.25 mag 
throughout, but we have checked that our conclusions also 
hold for AM = 0.5 or AM = 0.1. 

The bias of galaxies in the range M, ± AM/ 2 is then ap- 
proximated as the number-weighted average clustering of ha- 



los of mass irii ± Am/2. We can express bi as 



bi-- 



m,-+Am/2 



b(m)— — am 
am 





nnij+Am/l 


X 






J m/-Am/2 



dnh 



/ m—Amjl 

where b(ni) = b\v(rn)\ as defined in Equation[8] 



dm 



-Am 



(12) 



2.4. Sample Covariance from Galaxy Abundance and 
Clustering 

Equation |9] provides the average expected abundance of 
galaxies in a survey, given our fiducial luminosity function 
model ^(M). Owing to spatial density fluctuations on large 
scales, the actual measured number density of galaxies in the 
magnitude range M± AM/2 at location x will be 



n;(x,z) = n/[l-l-/7,(5(x,z)] 



(13) 



where (5(x,z) is the local linear overdensity and the bias bi is 
determined in Equation [121 The large scale structure of the 
matter density field will cause the galaxy counts to covary. 
The sample covariance Sij between galaxies in the /-th and 
/■-th magnitude bins is simply the average squared difference 
between the measured galaxy density n and the expected av- 
erage galaxy density h for each bin. We can then write the 
sample covariance as 



5,7 = {(ni-ni)(nj-nj)),. 



(14) 



where the average is take over all A^fieids fields of the survey. 
Given Equations [Tl [51 l9l [T2l and [131 we can evaluate the ele- 
ments of the sample covariance matrix S as 



fields 



(Pk 

(2^ 



Wiik)W*(k)P(k) (15) 



where W,(k) is the A:-space window function for the survey 
field volume of the i-th magnitude bin. Depending on, e.g., 
the redshift distribution of sources with different magnitudes, 
or some luminosity-dependent completeness, we could have 
W,(k) y W/k) in general. However, throughout the rest of the 
paper we will consider only galaxy densities and variances 
within the entire effective survey volume, such that the ele- 
ments of the sample covariance matrix S refer to luminosity 
bins within the same volume of each field. We will therefore 
write W;(k)W;(k) = |Wv(k)p, where 



Wv(k) = sine ( 1 sine I -y^ 1 sine I — 



(16) 



is an approximate A:-space window function for the effec- 
tive volume y of a survey at comoving radial distance r, 
comoving radial width 6r, and rectangular area 8,i x 8y in 
square radiansl3 The function sinc(x/2) = 2sin(ji:/2)/x is the 
Fourier transform of the Heavisid e Il(x) unit box. Similar 
win dow functions wer e adopted by Newman & DavisI ( I2002h 
and lStark et alJ (12007b) in their estimates of cosmic variance. 

Unless otherwise specified, when discussing the cosmic 
sample variance uncertainty or error we will refer to the aver- 
aged quantity 



CTCV = {b) D{z)auM / V A^fields , 



(17) 



* A rectangular survey will have a larger on-sky footprint at r + 5r than 
at r. For 0, and Sy of interest to this paper (< 1 deg) and 5r/r <g 1, we 
have checked that bo th W (x) and W(k) are well-approximated by the window 
function in Equation ll6l and its transform. 
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where {b) is the average bias of all galaxies in a survey field. 
This quantity ctcv is the cosmic variance uncertainty that is 
often reported for surveys, but is distinct from the elements 
sample covariance matrix Sij since the latter involves the bias 
of galaxies in individual luminosity bins. 

2.4. 1. Multiple Fields and Sample Covariance 

Splitting a survey into A^seids multiple fields can reduce the 
sample covariance in the combined data by probing statis- 
tically independent regions in space. Th e sample variance 
will scale roughly as S oc 1 /A^fieWs (e.g., iNewman & DavisI 
1^02); however, the actual gain depends on the shape of the 
ACDM power spectrum through the survey volume and ge- 
ometry. For a fixed amount of observing time, splitting a sur- 
vey into A^fieids = 2 fields will rescale the volume of each field 
by y cx 1 /A^fieids and result in a corresponding increase in the 
typical dark matter density fluctuations in each field. For very 
large surveys, the decrease in the volume per field can (at least 
partially) offset the gains achieved by p robing multiple inde- 
pendent samples (see also Munoz et al. 2009). 

The left panel of Figu re [T] shows the RMS density fluctua- 
tions D(z)(7dm / V-^fieids in a survey at redshifts 6.5 < z < 7.5 
as a function of total area for multiple fields (A^fieids = 1,2,4). 
For a galaxy survey, the sample variance in each luminos- 
ity bin will be increased by a factor of the galaxy bias (see 
Equation [TSl l. While the uncertainty from RMS density 
fluctuations will improve with the addition of statistically- 
independent samples, the fractional improvement is less than 
1 - l/V-^fieids for large volumes. The right panel of Fig- 
ure [T] shows the fractional improvement gained by multiple 
fields. For a flat power spectrum, the improvement would 
be 1 - \/y/2 = 0.293 for A^fieids = 2 and 1 - 1 /V4 = 0.5 for 

A^fields = 4. 

2.5. Poisson Variance from Galaxy Abundance 

Number-counting statistics will naturally introduce a Pois- 
son variance into the galaxy number count statistics. The di- 
agonal Poisson variance matrix P will only add to the total 
covariance for counts within a single magnitude bin (i.e., only 
when / = j). For definiteness, we will express the Poisson co- 
variance as 



V, 



(18) 



where the Kronecker 5ij = 1 for / = j and 5ij = for / ^ j. 

3. PARAMETER ESTIMATION AND THE FISHER MATRIX 



iFisheij (Il935h illustrated how to infer inductively the prop- 
erties of statistical populations from data samples. By approx- 
imating the likelihood function as a Gaussian near its maxi- 
mum and assuming a parameterized model, the uncertainties 
in the model parameters allowed by a future data set can be 
estimated directly from the data covariances. Interested read- 
ers should refer to the excellent and detaile d discussion of 
the Fisher matrix formalism provided in §2 of 'T egmark et alj 
([T997), but we outline the general approach below. 

We aim to estimate the uncertainties on model parameters 
(the "parameter covariance matrix" C) achieved by the 
data produced by some fiducial survey. The quality of the 
future data for each fiducial survey will be characterized by 
the "data covariance matrix" D. The elements Dij of the total 
data covariance matrix are simply the sum of t he sa mpl e co- 
variance and Poisson uncertainties described in ^2.4l and ^2.51 



which we can write as 

D;j = Sij + Pij, (19) 

where the P, , only contribute when / = j (see Equation [TS]). 
Following Lima & Hd (120041) . who applied the Fisher matrix 
approach to the parameter estimation of the mass-observable 
relation in galaxy cluster surveys, we will express our approx- 
imate Fisher matrix as 



dni 



D 



V dtii 1 



D" 



as 



-D 



dS 



dPtj. dpy 



(20) 



(see, e.g., iHold er etalJ I2OOII: iHu & Kra vtsovl I200I 
Lima&Hu 2005: iHu & CohnI |2006|; ICunha&Evri^ 
2009; Wu et al. 2009, and especially the discussion in §111 of 
Lima & Hu 2004). Here, Tr(A) = Yl7=i^ii f^'" ™ ^ nia- 
trix. The vector elements p^ reflect the parameters of the data 
model n. The first term of Equation|20]models second deriva- 
tives of t he likelihood function in the Poisson error-dominated 
regime (iHolder et al.l 1200 1 ), while the second term models 
the s ample covariance-dominated regime (Tegmark et al] 
119971 see also Appendix A of Vogeley & Szalay 1996). 
The derivatives dn/dp/^i of the luminosity function model 
are computed directly by differentiating Equation |9] The 
derivatives dS/dp^ of the sample covariance matrix are 
evaluated numerically since changes to the model luminosity 
function alter the galaxy bias b in Equation [TSl for a given 
luminosity bin in a nontrivial way. 

Once the Fisher matrix F is calculated, estimating the pa- 
rameter covariance matrix C becomes straightforward. The 
elements of the parameter covariance matrix are approxi- 
mated as 



The marginalized uncertainty on parameter is then 



'2 



(21) 



(22) 



Similarly, we can estimate the unmarginalized error on each 

-I /2 

parameter as a'^ = Fi_i,_i . However, in what follows when we 
discuss the "error" or "uncertainty" on luminosity function 
parameters we mean the marginalized error unless otherwise 
stated. 

We will apply the above formalism to estimate the relative 
constraining power of possible galaxy surveys, but we will 
focus on evaluating such surveys in the context of existing 
and forthcoming data from observational programs already 
underway (i.e., the WFC3 UDF GO and GOODS ERS data). 
Our statistical formalism provides a simple way to incorporate 
constraints from prior data. The combined constraints Ccombo 
of a prior observation Cprior supplemented by the forecasted 
constraints of a future survey C can be estimated as 



-^combo ■ 



(c-'+c 



priory 



(23) 



or, in other words, the combined parameter covariance matrix 
is the inverse of the sum of the Fisher matrices of the prior 
and future surveys. Depending on the magnitude of the off 
diagonal elements of C and Cpiim, the combined parameter 
covariance matrix Ccombo can provide a substantially different 
correlation between parameters than either the prior or future 
surveys produce individually. As a result, the marginalized 
uncertainty on parameters can benefit substantially by com- 
bining surveys with different characteristics. These ramifica- 
tions of Equation|23]will become more apparent in ^ 
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Fig. 1 . — Root-mean-squared (RMS) density fluctuations in a survey at redshift 6.5 < z < 7.5 as a function of total area (right panel). Uncertainties owing to 
RMS density fluctuations scale with the product of the growth function D(z) and the z = RMS dark matter fluctuations ctdm- If the ACDM power spectrum 
""DM for the rectangular survey geometry was flat, multiple field surveys would improve their combined RMS density fluctuations by 1 / VMiclds if they probed 
widely-separated, statistically-independent regions. However, the shape of the ACDM power spectrum ctdm limits this improvement since the multiple fields 
each probe a volume V /Nums and the power increases toward small scales. Shown are the RMS density fluctuations for multiple field surveys (Wflcids = 1, solid 
line; A'acMs = 2, dashed line; Wfieids = 4, dotted line). The left panel shows the fractional improvement in uncertainties owing to density fluctuations gained by 
splitting survey into Wfleids = 2 (dashed fine) or Wjeids = 4 smaller fields of equivalent total area. If ctdm were independent of scale, the fractional improvement 
would be a constant 1 - 1 / ^2 = 0.293 for Wfleids = 2 and 1 - 1/^4 = 0.5 for Nf^M, = 4. 



4. MODEL FOR SURVEY DATA 

Our Statistical formalism for forecasting constraints on 
properties of the galaxy population requires a model for the 
survey data. Given the approach outlined in ^ the relevant 
characteristics of each survey include the total area AjQt, num- 
ber of fields A^fieids, and the minimum and maximum redshifts 
of the survey Zmin and z^ax (that, in combination with A^ot, de- 
termine the survey volume V). The limiting magnitude depth 
of the survey Mmax strongly influences source statistics of the 
survey by determining the faintest luminosity bin calculated 
via Equation |9] The completeness of the survey /comp and 
halo occupation fraction /occ change the cosmic sample vari- 
ance by altering the halo-galaxy correspondence in Equation 

m 

Of these survey characteristics, we will keep Zmin, Zmax, 
/comp, and /occ fixed between surveys. We will assume that 
the surveys are effectively volume-limited (/comp =1) over 
the redshift range of interest. Given a complete volume- 
limited survey, the choice of minimum and maximum red- 
shifts roughly corresponds to the filter choice defining a 
dropout selection. We will adopt Zmin = 6.5 and Zm^x = 7.5, 
which roughly approximates the redshift sele ction of the (zhsq- 
Yios) vs. (F105--/125) color selection of Oes ch et all (l2010bt 
see their Fig. 1). Similar selections can be defined for /814- 
dropouts. Our calculations can be easily extended to differ- 
ent redshift selections, but we adopt this redshift range since 
the fiducial abundance of z ~ 7 WFC3 UDF GO galaxy can- 
didates app ears increasingly r obust (see, e.g., the discussion 
in §2 of McLu re et alj 1201 Oh . the characteristic ultraviolet 
(UV) luminosit y of galaxies is decreasing with redshift (e.g., 
iBouwens et al.l l2008), and the abundance of dark matter halos 
hosting galaxies is rapidly declining at earlier epochs. 

Our model surveys will consist of WFC3 //-band cover- 
age with equal coverage in an additional, bluer WFC3 filter 
The existing and ongoing UDF GO and GOODS ERS sur- 
veys will use the F160W band (see ^below), but using the 



F14QW band buys w 0.3-0.5 magnitudes in sensitivity for 
the same exposure time, depending on the source luminosity 
(see below). While the dropout color selection is perhaps bet- 
ter for F160W, we will assume in our forecasts that future 
surveys will utilize F14QW. Our results would be similar for 
F16QW surveys to similar limiting depths. We will charac- 
terize the abundance of high-redshift galaxies in terms of a 
the rest-frame ultraviolet (UV) luminosity function. We must 
therefore adopt a color conversion between //ab magnitude 
and rest-frame luminosity appropria te for z ~ 7. I n an app rox- 
imation to the conversion used by lOesch et alJ (l2010bh . we 
estimate that //ab ~ 29.0 translates to Muv ~ -18.2. We have 
checked that our general conclusions about the relative con- 
straining power of survey designs are insensitive to changes 
in this conversion (e.g., ±0.4 magnitudes in Muv)- 

Lastly, HST observations are conducted using some num- 
ber A^oibits per pointing that effectively determines M^^^. Dur- 
ing each orbit the field visibility depends on the field decli- 
nation, and the available on-source integration time also de- 
pends on observatory and instrument overheads such as guide 
star acquisition, filter changes, dithering, and readout. The 
UDF GO and GOODS ERS surveys aie at a declination of 
S Ri -27 deg, and for ease of comparison we will assume all 
future surveys have \S\ < 30 deg. This roughly equatorial 
declination range provides a visibility of 54 minutes/orbifl 
Given the additional observatory and instrument overheads, 
we will calculate all sensitivities using a 46 minutes/orbit ex- 
posure time. Given the compact ch aracter of the obse rved 
WFC3 z ^ 7 galaxy candidates (e.g.. lOesch et al.|[2010a ^. we 
will report optimum 5-cr point source sensitivities. Table [T] 
lists these sensitivities for a flat /v spectrum source, as a func- 
tion of A^oibits for both the FUOW and F160W filterfl 



^ see Table 6.1 of|http://www.stsci.edu/hst/proposing/docunients/piimer/| 
^ Computed using the WFC3 IR Channel Exposure Time Calculator, 
http://etc.stsci.edu/webetc/mainPages/wfc31RlmagingETC.jsp 
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TABLE 1 

Optimum 5-ct Point Source Sensitivity 
vs. Exposure Time Per Pointing 



A'oibits 


//i4o [AB Mag.] 


//,60 [AB Mag.] 


0.5 


27.00 


26.62 


1.0 


27.43 


27.07 


2.0 


27.84 


27.49 


3.0 


28.07 


27.73 


4.0 


28.23 


27.89 


6.0 


28.46 


28.12 


8.0 


28.62 


28.28 


19.0 


29.10 


28.76 


38.0 


29.47 


29.14 


125.0 


30.12 


29.78 



5. EXISTING SURVEYS 

The discussion in ^ makes clear that the combination 
of different survey designs can potentially provide increased 
constraints beyond that achieved by individual data sets. Even 
duplicate surveys will reduce the Poisson variance and poten- 
tially the sample variance (especially if the fields are widely 
separated on the sky). In the absence of significant system- 
atic biases, using prior data will generally improve the ex- 
pected parameter uncertainty obtained by future experiments. 
We will therefore rely on the expected constraints achieved 
by the ongoing WFC3 UDF GO (PI Illingworth, Program ID 
1 1563) and GOODS ERS (PI O'Connell, Program ID 1 1359) 
programs to augment the fiducial survey designs evaluated in 
^and ^ In this section, we will calculate the expected con- 
straints provided by the UDF GO and GOODS ERS data. 

Table|2]describes the field geometry, number of fields A^fieids, 
total area, and expected //-band 5-0- point source depth for 
the UDF GO and GOODS ERS survey designs. Numerous 
analyses of the initial UDF GO data release have already been 
performed (e. g., Bouwens et al. 2010a; Oesch et al. 2010b3; 
McLure et al. 2010; Bunker et al. 2009; Jvlc Lure et all 120101: 
Yan et al. 2009; Finkelstein et al. 2009), but we will consider 
the expected constraints provided by the entire 192 orbit pro- 
gram. The GOODS ERS data has not yet been released (but 
for initial analyses on u nreleased ERS data see IWilkins et alj 
IMotlLabbe etsIM o§). and we will also use the Fisher ma- 
trix approach to estimate the constraints provided by that sur- 
vey. 

The UDF GO program is comprised of three WFC3 point- 
ings. Two of the WFC3 pointings use the F160W filter with 
19 orbits. Using the WFC3 ETC, we estimate these obser- 
vations will reach Hab ~ 28.76. The UDF GO observations 
also will have 38 F160W orbits in the HUDF that will reach 
HAB~29.l4^ The remaining 116 orbits in the program will 
be used for observing in bluer filters. 

The GOODS ERS survey will have 3-orbit depth in F160W, 
using 24 orbits (out of a total 104) for //-band observations. 
We estimate that these observations will reach a sensitivity of 
//ab ~ 27.7^ The remaining 80 orbits in the program will be 
used for observations with other filters and grisms. 

5.1. Forecasted Constraints for Existing Surveys 



^ The UDF GO proposal e.stimates their sensitivitie.s as Hab ^ 28.6 for 19 orbits and 
//ab ^ 29 for 38 orbits. See http://www.stsci.edU/observing/phase2-public/I I563.pdf 

^ The GOODS ERS proposal estimates the 3-orbit sensitiv- 
ity using 40 minutes/orbit integration as //ab ^ 26.5. See 
lhttp://ft^w.stsc i.edu/hstyproposing/old-proposing-files/goods-cdfs.pdf and 
"http://www.stsci.edu/obsei'ving/phase2-public/l I359.pdf Our method for estimating 
the sensitivity would provide //ab ^ 27.65 for 40 minute,s/orbit. 



The forecasted constraints achieved by the UDF GO and 
GOODS ERS surveys are plotted in Figure |2] Each panel 
shows the constraints for the UDF GO (blue region) and 
GOODS ERS (light blue region) surveys separately, and in 
combination (dark blue region). The constraints are plotted 
for the (left panel), M*-a (middle panel), and 0+-a 

(right panel) projections. In Figure |2] (and in similar figures 
throughout the paper), the shaded regions correspond to stan- 
dard Gaussian contours. 

Figure |2] highlights some general properties of the perfor- 
mance of different kinds of surveys for providing luminosity 
function parameter constraints, as well as specific features of 
the UDF GO and GOODS ERS surveys: 

• The covariances between luminosity function parame- 
ters are significant and positive. In terms of the Pear- 
son's correlation coefficient p for parameters x and y, 
defined as 



C„ 



(24) 



the forecasts calculate typical correlation coefficients of 
p > 0.9. For a given a narrow range of logjQ i/i* or 
a are permitted by the data, even as the marginalized 
uncertainties can be as large as ~ 30-40% fractionally 
for log J and a. 

• The orientation of the constraint ellipse forecasted for 
each survey can differ significantly depending on the 
parameter uncertainties, even if the parameter corre- 
lation coefficients for the separate surveys are similar 
The orientation of the constraint ellipse major axis with 
respect to the x-axis in the x-y parameter plane can be 
characterized by the angle 



9 = - arctan 
2 



IpUxCfy 



1 

- arctan 
2 



2C,. 



C -C 



(25) 



which depends on the correlation p and the parameter 
uncertainties. If the angle Q differs between separate 
surveys, then the constraints achieved by combining the 
surveys can improve dramatically. 

For reference, the calculated marginalized and unmarginal- 
ized errors for p as well as the Pearson's correlation coeffi- 
cient p and the angle O for each pair of parameters are listed 
for the existing surveys in Table|4] 

As Figure |2] and Table H show, the UDF GO and 
GOODS ERS surveys will akeady provide interesting con- 
straints on the abundance of z ~' 7 galaxies. When combined, 
the unmarginalized uncertainties on the LF parameters will 
be AM* - 0.1 mag, Alogm^* < 0.1, and Aa 0.15. For 
the full UDF GO survey, we find that the ummarginalized un- 
certainty for the faint-end slope is Aa ^ 0.16. Using a sin- 
gle UDF GO field an d a limiting depth of F16QW - 29 AB 
for a single pointing, lOesch etal] (120 1 Obi) report a faint-end 
slope uncertainty of Aa ~ 0.33 when and M* are fixed 
(i.e., the unmarginalized uncertainty on a). If we use the 
same single pointing area and depth, and the same cosmol- 
ogy, our estimate of the unmarginalized uncertainty would in- 
crease to Aa - 0.26. The GOODS ERS and UDF GO sur- 
veys are complementary in that the depth of the UDF GO sur- 
vey provides a beneficial constraint on the faint-end slope a. 
This UDF GO constraint on a rotates the UDF GO error el- 
lipse relative to the GOODS ERS constraint in the M* - a 



8 



TABLE 2 
Existing Surveys 



Survey 


Field Geometry 


A'ficlds 


Total Area 


Worbits" 


//-band Depth 


Ref. 


Name 


[Arcmin. X Arcmin.] 


[#] 


[Sq. Arcmin.] 


[#] 


[AB Mag.] 




UDFGO 


2.05' X 2.27' 


2 


9.3 


19 


28.76 


1 




2.05' X 2.27' 


1 


4.7 


38 


29.14 




GOODS ERS 


5.2'x iQ.3"' 


1 


53.3 


3 


27.73" 


2 



References. — (1) http://www.stsci.edu/observing/ phase2-public/11563.pdf (2) 
http://www.stsci.edu/hst/proposing/old-proposiiig^fiIes7g6ods-cdfs.p3i: 
■' Number of orbits per pointing. 

^ The GOODS ERS survey is a 2 X 4 WFC3 IR mosaic dithered to match the Ultraviolet/Visible 
channel field of view. 
For details, see the discussion in ^ 
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Fig. 2. — Forecasted constraints o n z ~ 7 lu minosity function parameters expected from the forthcoming UDF GO and GOODS ERS WFC3 data. The 
constraints are calculated for a Schechte3 U976I) luminosity function with a chai'acteristic luminosity Mj,, normalization (f)^,, and faint-end slope a. Shown are 
the l-cr constraints in the M* — logjQ 0* (left panel). Mi, -a (middle panel), and logjQ (p-t — a (right panel) space projections for the UDF GO (blue region) and 
GOODS ERS (light blue region) surveys. Also shown are the constraints expected by combining both surveys (dark blue region). The depth of the UDF GO survey 
will provide a better constraint on the faint-end slope than GOODS ERS, but the differences between their parameter covariances make them complementary. 



and 0* - a projections, thereby reducing the corresponding 
parameter uncertainties. Individually, the GOODS ERS un- 
certainties will be considerably larger than those obtained by 
the UDF GO survey, since the GOODS ERS survey lacks suf- 
ficient depth to tightly constrain the LF faint-end slope and is 
not wide enough to tightly constrain or 

While the UDF GO and GOODS ERS surveys achieve 
appreciable unmarginalized constraints, the covariances be- 
tween the LF parameters are large. The marginalized uncer- 
tainties calculated for the LF parameters are AM,^ ^ 0.5 mag, 
A(l>^ ^ 0.4, and Aa ^ 0.4. Accounting for covariances these 
marginalized parameter uncertainties correspond to a frac- 
tional uncertainty in the total number of galaxies with Muv < 
-18 of « 25%, increasing to a factor of w 2 uncertainty in 
the total number of galaxies with Muv < M^,. To improve 
the constraints on the number of galaxies with Muv < -18 
(Muv < M*) to « 5% (~ 50%) would require marginalized 
parameter uncertainties of approximately AM,^ ^ 0.2 mag, 
A(pi, ~ 0.2, and Aa ^ 0.2 depending on their covariances. To 
reach such constraints, these UDF GO and GOODS ERS sur- 
veys would need to be complemented by either wider area or 
deeper surveys. We now consider some fiducial model sur- 
veys that could achieve these constraints in combination with 
the UDF GO and GOODS ERS data. 

6. MODEL SURVEYS 

The complete UDF GO and GOODS ERS surveys will pro- 
vide extremely interesting initial data on the abundance of 
z ~ 7 galaxies, but the marginalized uncertainties on the LF 



parameters achieved by those surveys will still permit uncer- 
tainties of ^25% in the total number of galaxies at Muv S 
-18. We can repeat the calculations from ^ for fiducial 
model surveys to illustrate what constraints wider or deeper 
surveys can achieve when combined with the UDF GO and 
GOODS ERS data. 

The model surveys are designed to be appropriate for a HST 
Multi-Cycle Treasury Progranfl which can receive up to 750 
orbits per HST cyclq3 We consider six possible model sur- 
veys that we estimate would require 450 - 900 total orbits to 
acquire filter coverage with the WFC IR channel. 

As discussed in ^ we will assume the model surveys will 
use the F140W filter owing to its increased throughput rela- 
tive to F16QW. The sensitivity of each survey is determined 
by selecting a number A^orbits of orbits per pointing, assuming 
46 minutes/orbit exposure time, and using the WFC3 IR chan- 
nel ETC. The total number of orbits for each survey were then 
determined by selecting the number of fields A^fieids, a mosaic 
geometry per field, and multiplying the number of pointings 
in each mosaic by A^orbits (and then doubling to account for 
comparable coverage in a bluer WFC3 filter). 

The survey models are designed to cover a large range in 
total area (Atot ~ 14-3600 square arcmin), field numbers 
(A^fieids = 1-4), orbits per pointing (A^orbits = 0.5- 125), limiting 
depth (Hab ~ 27-30), and total number of orbits (450-900). 
We design each survey to approximate possible HST WFC3 

' See, e.g., http://www.stsci.edu/institute/org/spd/mctp.html/| 

See http://www.stsci.edu/institute/org/spd/HST-multi-c ycle-treasury1 
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tilings of existing surveys; as such, these model surveys repre- 
sent realistic extensions of existing HST and Spitzer surveys 
to hundreds of orbits of WFC3 coverage. Summaries of the 
model surveys can be found in Table |3] and are ordered by 
decreasing limiting depth and increasing total area. Brief de- 
scriptions of the models follow: 

6.1. Survey A 

The performance of the UDF GO survey suggests that an 
interesting possible survey design would be a set of narrow 
pencil beam surveys with sufficient depth to reach a few nJy 
sensitivity. Extending each of the three single WFC3 point- 
ing UDF GO fields to - 125 orbits in F140W would achieve 
Hab ~ 3O.l03- For surveys of ~ 10 square arcmin total area, 
using A^eids > 1 results in a substantial reduction of sample 
variance (see Figure [T]i- The model Survey A will therefore 
use A^fieids = 3, A'oibits = 125, and Atot = 13.96 square arcmin 
(three WFC3 pointings), and 750 total orbits including cover- 
age in a bluer WFC3 filter. For calculating constraints from 
a combination of Survey A with existing data, we assume the 
Survey A fields will be able to leverage the GOODS ERS data 
but will duplicate the UDF GO datO 

6.2. Survey Bl 

Another template for a model survey is deep WFC3 cover- 
age of the GOODS survey fields. A 5 x 7 WFC3 mosaic could 
cover a field of size 10.3 ' x 15.9', similar to the GOODS fields 
dGiavalisco et al.l 120041) . Covering A^tieids = 2 fields the size 
of the GOODS fields would require 70 pointings, and would 
cover a total area of Atot = 326 square arcmin. Using A^oi-bits = 8 
orbits per pointing would reach //ab = 28.6 in F140W, and 
would require a total of 896 orbits (including coverage in a 
bluer WFC3 filter). Survey Bl is the most expensive survey 
we consider. For calculating combined constraints utilizing 
existing, we will combine Survey B 1 with the UDF GO data 
but ignore the duplicated GOODS ERS F160Wdat£O 

6.3. Survey B2 

To gain intuition about the relative value of depth and area 
for constraining high-redshift galaxy populations, we will 
consider variations of the GOODS-like survey. Survey B2 is 
identical to Survey B 1 in number of fields and pointings, but 
would achieve a reduced depth of A^orbits = 6 orbits per point- 
ing (//ab = 28.46). The total number of orbits required for 
Survey B2 is 672 (including equal coverage in a bluer WFC3 
filter). When determining combined constraints with existing 
data, we will combine Survey B2 with the UDF GO survey. 

6.4. Survey B 3 

Same as Survey Bl and Survey B2, but to A^oi-bits = 4 orbits 
per pointing (//ab = 28.2) depth. Survey B3 would require 
448 total orbits. For combined constraints with existing data, 
we will combine Survey B3 with UDF GO. 

' ' Surveys this deep require many WFC3 frame exposures to avoid image 
persistence. We ignore the impact of any additional related overhead on the 
available exposure time. 

The F160W data from UDF GO could be incorporated the same //ab- 
band depth, which could potentially decrease the total orbits for this survey 
by ~ 150. Our general conclusions are not strongly influenced by choosing 
this alternative. 

" Most of the additional constraint achieved by combining with existing 
data comes from the UDF GO data, so this choice is not critical for our gen- 
eral conclusions. However, using the GOODS ERS data could potentially 
decrease the required orbits for a GOODS-like survey by ~ 24 orbits. 



6.5. Survey C 

An existing survey with a combination of large area and 
infrared depth is the Spitzer Extended Deep Survey (SEDS; 
iFazio et al.ll2008b . which was designed to cover 0.9 deg over 
five fields to 12 hour/pointing depth with the warm Spitzer 
Infrared Array Camera 3.6/im and 4.5/im channels. Exactly 
reproducing the SEDS survey with WFC3 would likely be 
prohibitively expensive, so we will instead consider a feasi- 
ble WFC3 survey with a design similar in spirit to SEDS. 
Our SEDS-like Survey C will consist of A^fieids = 4 fields of 
4x13 pointing mosaics (each of size 8.2' x 29.5'), for a total 
area Atot = 967.6 square arcmin. A depth of A^oi-bits = 2 orbits 
per pointing (//ab = 27.8) would then require 832 orbits (in- 
cluding equal coverage in a bluer WFC3 filter). For calculat- 
ing combined constraints with existing data, we will combine 
Survey C with both the UDF GO and GOODS ERS fields. 

6.6. Survey D 

The largest HST survey to date is th e equato rial Cosmic 
Origins Survey (COSMOS Scoville et al l2007bl) . which cov- 
ers 2 deg^ with the ACS /-band. As with the SEDS-like 
Survey C, exactly reproducing the COSMOS survey with 
WFC3 would likely be prohibitively expensive. Instead, we 
consider a 1 deg^ (Atot = 3629 square arcmin) survey with a 
single 26 x 30 mosaic (59.0' x 61.5') to A^orbits = 0.5 orbits 
per pointing (//ab = 27) depth. Survey D is the widest and 
shallowest design we consider, and would require 780 orbits 
to complete (including equal coverage in a bluer WFC3 fil- 
ter). We will combine Survey D with both the UDF GO and 
GOODS ERS surveys for purposes of calculating combined 
constraints incorporating existing data. 

6.7. Field Size Comparison 

We show an illustrative comparison of the existing and 
model survey areas in Figure|3] The UDF GO, GOODS ERS, 
and Surveys A, Bl, B2, B3, C, and D areas are shown as white 
boxes overlaid on a thin 10/;~'Mpc slice through a ACDM 
cosmological simulation of comoving size L = 250/!~'Mpc at 
z ^ 7 (Tinker et al. 2008). The blue scale image shows the 
projected dark matter surface density calculated from the dark 
matter particle distribution of the simulation. The comoving 
length scale corresponding to an angle of 9 = \ deg at z = 7 
is 108.7/i"'Mpc for the adopted WMAP5 cosmology. This 
comparison illustrates the characteristic angular size of large 
scale structures at z ^ 7, as well as the survey areas required to 
probe representative samples of the high-redshift dark matter 
density distribution. The separation between fields is not to 
scale, and model surveys incorporating different fields would 
likely be more widely spaced to probe statistically indepen- 
dent regions on the sky. 

7. FORECASTED CONSTRAINTS FOR MODEL SURVEYS 

The forecasted constraints calculated for the model Sur- 
veys A, Bl, B2, B3, C, and D are summarized in Table |5] 
and presented in Figures |4]|9] In each figure, the shaded ar- 
eas show the projected constraints for each model survey in 
the Mi, - (t>i, (left panel). Mi, - a (middle panel), and 
(right panel) LF parameter planes. The axes ranges in Figures 
|4l|9]are identical (and much smaller than in Figure|2]i, and the 
plotted constraints are directly comparable. A description of 
the forecasted constraints for each model survey follows; 
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TABLE 3 
Model Surveys 



Model 


Field Mosaic 


Field Geometry 




Total Area 




//-band Deptli 


Total Orbits'' 


Name 


[# Point. X # Point.] 


[Arcmin. X Arcmin.] 


[#] 


[Sq. Arcmin.] 


[#] 


[AB Mag.] 


[#] 


Survey A 


1 X 1 


2.05' X 2.27' 


3 


13.96 


125 


30.12 


750 


Survey B 1 


5x7 


10.3'x 15.9' 


2 


325.7 


8 


28.62 


896 


Survey B2 


5x7 


10.3'x 15.9' 


2 


325.7 


6 


28.46 


672 


Survey B3 


5x7 


10.3'x 15.9' 


2 


325.7 


4 


28.23 


448 


Survey C 


4x 13 


8.2'x 29.5' 


4 


967.6 


2 


27.84 


832 


Survey D 


26 X 30 


59.0' X 61.5' 


1 


3628.5 


0.5 


27.00 


780 



^ Number of orbits per pointing. 

We assume each survey will require comparable coverage in two WFC3 filters, which doubles the required number of total orbits. 



ACDM, z~7 f. 



. □ □ □ UDF/Surv( 
,* 



Survey G ' 




: 250 Comoving h" Mpc, d = 2.3 degrees at z= 



Fig. 3. — Existing and model survey areas compared with the large scale dark matter structure at z ~ 7. Shown are the UDF GO, GOODS ERS, Sur- 
veys A, Bl, B2, B3, C, and D areas (white boxes), projected onto a surface density map of a thin 10/i^'Mpc slice through a ACDM cosmological simulation 
of size L = 250/)"' Mpc ( Tinker et alj i2008y The number of fields and survey areas of UDF GO and Survey A are identical. This comparison illustrates the 
characteristic angular size of large scale structures at z ^ 7, as well as the survey areas required to probe representative samples of the high-redshift dark matter 
density distiibution. The separation between fields is not to scale, and model surveys incorporating different fields would likely be more widely spaced to probe 
statistically independent regions on the sky. 



7.1. Survey A 

Figure |4] shows the forecasted constraints for Sur- 
vey A (Atot = 14 square arcmin, A^fieids = 3, Hab = 
30.1, blue region), GOODS ERS (Hght blue region), and 
GOODS ERS and Survey A combined (dark blue region). 
Survey A could find > 250 z ^ 7 galaxies to H^b ^ 30.1, 
with a Poisson variance in the galaxy count of w 6%.. The av- 
erage galaxy bias for this depth and area is (b) w 5.3, which 
results in a sample cosmic variance of ctcv ~ 0.16 fraction- 



all\Pl The constraints achieved by Survey A are therefore 
cosmic variance limited. Survey A is the deepest model sur- 
vey we consider, and results in the tightest forecasted con- 
straints on the LF faint-end slope (Aa « 0.1, marginalized). 
GOODS ERS complements Survey A by providing an im- 

'"* For our cosmology and the lOesch et al J )2010bl) estimated LF, the co- 
moving volume of a single WFC3 pointing to //ab ~ 29 depth at the red- 
shifts of interest would have a cosmic variance uncertainty of ctcv ~ 0.32 
(c.f . Oesch et al. 2010b). 
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proved combined constraint on the LF normalization (A(/)* w 
0.15) and characteristic magnitude (AM* w 0.22). The com- 
bined Survey A and the GOODS ERS survey also produce 
relatively low correlation coefficients (p w 0.6-0.85) com- 
pared with other surveys combinations; the orientation of 
the - a joint constraint from Survey A is only slightly 
inclined (6 = 12 deg), and allows the GOODS ERS survey 
(6 = 55 deg) to improve the combined constraints on M*. 

7.2. Surveys Bl, B2, and B3 

These model surveys illustrate how increasing the survey 
limiting depth over a moderate area alter the forecasted LF pa- 
rameter constraints. These surveys share a common total area 
(Atot = 326 square arcmin) and number of fields (A^fieids = 2), 
but span a factor of two in integration time (sa V2 in sensitiv- 
ity, //ab = 28.2-28.6, see Table |3]l. Because the extra depth 
probes more abundant, lower-luminosity galaxies, the typi- 
cal galaxy bias ({b) « 6.5) and cosmic variance uncertainty 
(ccv ~ 0.142) in Survey Bl would be smaller than for either 
Survey B2 ({b) « 6.7, ctcv ~ 0.146) or Survey B3 ((/?) « 7.0, 
(Tcv ~ 0.152). Similarly, the extra depth affords more ob- 
served galaxies (N ^ 1 100) and less Poisson uncertainty (3% 
fractionally) for Survey Bl than for Survey B2 (A^ w 850, 
3.4%) or Survey B3 (A^ w 610, 4.1%). 

The Fisher matrix calculations translate the Poisson and 
cosmic variance uncertainties into constraints on the LF pa- 
rameters, and Figures|5]|7]show how the parameter constraints 
scale with limiting magnitude for Surveys Bl, B2, and B3. In 
each figure, the shaded areas show the constraints achieved 
by UDF GO (light blue), the model surveys (blue), and the 
combination of UDF GO and each model survey (dark blue). 
The LF parameter constraints are also listed in Table |5] 

Of these three surveys. Survey Bl achieves the best 
combined parameters constraints (AM* = 0.21, Alog^^* = 
0.176, Aa = 0.172). However, the relative gain over Sur- 
vey B2 (AM* = 0.24, Alogio0* = 0.19, Aa = 0.21) and Sur- 
vey B3 (AM* = 0.25, Alogioc/)* = 0.21, Aa = 0.25) are rel- 
atively modest (20% improvement in AM* and Alogjoc/)*, 
and 40% in Aa). Most of the relative improvement owes 
to the increased constraint on the LF faint-end slope for Sur- 
vey Bl, since the three surveys are essentially identical for 
galaxies with //ab < 28.2, have a similar orientation of their 
error ellipse in the M*-0* projections (9 w 39 deg), and have 
similar correlations between LF parameters. Combining with 
UDF GO results in a larger relative improvement in the LF pa- 
rameter constraints for Survey B2 (10%) and Survey B3 (20- 
25%) than for Survey Bl (5%). 

7.3. Survey C 

The next widest model survey design is Survey C, with a to- 
tal area of Atot = 967.6 square arcmin to //ab = 27.8 depth over 
A^fieids = 4 fields. Such a survey would find « 940 galaxies 
at z 7, with an average bias of (b) = 7.5, cosmic variance 
uncertainty of tJcv ~ 0.1, and Poisson uncertainty of 3%. 

Figure [8] shows the constraints for Survey C (blue region), 
the combination of UDF GO and GOODS ERS (light blue 
region), and the combination of all three surveys (dark blue 
region). The larger area of Survey C allows for better or com- 
parable combined constraints on the LF characteristic mag- 
nitude (AM* w 0.19) and normaHzation (Alogjyc/)* w 0.15) 
than deeper surveys over smaller areas. Owing to its weaker 
constraint on the faint-end slope, the error ellipses provided 
by Survey C are more highly inclined in the M*-a (0 ~ 



50 deg) and 0*-a (Q w 59 deg) projections than the UDF GO- 
GOODS ERS combined constraints (38 and 44 deg). When 
combined with UDF GO and GOODS ERS surveys. Sur- 
vey C would provide among the tightest constraints of the sur- 
veys we consider (with Survey A providing better combined 
constraints on a and Survey D providing better constraints on 
M*, </i*, and a). 

Of additional interest for a design like Survey C is some 
measure of the benefit of having A^fieids = 4 for constraining 
the z ^ 7 LF compared to a single contiguous field. We note 
that changing Survey C to a single field of the same total 
area and aspect ratio results in essentially no change to the 
constraints on the LF parameters (a fractional change of less 
than 1%). The cosmic variance uncertainty does improve by 
^ 25% (see Figure[T]i from <tcv ~ 0.13 when increasing the 
number of fields from A^fieids = 1 to A^seids = 4, but this improve- 
ment has little net effect on the LF parameter constraints. The 
marginalized constraints on the LF parameters are sensitive to 
the Poisson errors of individual magnitude bins on the bright 
end of the LF, and the Poisson error is independent of A^eids 
for surveys of fixed total area. For magnitude bins that are 
Poisson-uncertainty dominated, the improvement in the cos- 
mic sample variance gained by increasing A^fieids therefore may 
not strongly influence end constraints on the LF parameters. 

7.4. Survey D 

The widest and shallowest model survey design consid- 
ered is the single-field Survey D (Atot ~ 1 deg^, //ab = 27, 
Mieids = !)■ This model survey would find w 570 galax- 
ies at z ~ 7, probing only galaxies brighter than M* with an 
average bias of (b) ~ 9 with a cosmic variance uncertainty 
of (Tcv « 0.11 (dominating over the Poisson uncertainty of 
4.2%). Figure|9]shows the constraints that would be achieved 
by the combination of UDF GO and GOODS ERS (light blue 
region). Survey D individually (blue region), and the combi- 
nation of all three surveys (dark blue region). 

The constraints achievable by Survey D individually 
are comparable to the constraints provided by combining 
UDF GO and GOODS ERS, but would require roughly three 
times as much additional telescope time to complete. How- 
ever, the combination of Survey D with both UDF GO and 
GOODS ERS produces the strongest joint constraint of any 
survey design we considered (AM* « 0.136, AlogjQi/)* « 
0.11, Aa w 0.20). The orientation of constraint provided 
by Survey D individually is inclined (9 = [29.5,65.4,76.2]) 
relative to the UDF GO-GOODS ERS combination (6 = 
[38.9,37.6,44.3]), and results in a relatively low correlation 
between the LF normalization and faint-end slope (p « 0.77). 
While other survey designs produce better constraints on the 
faint-end slope, the joint constraint region shown in Figure |9] 
produces an uncertainty in the LF that is better than w 6% at 
all relatively bright (//ab ^ 28) magnitudes. 

8. DISCUSSION 

We have considered the problem of forecasting constraints 
on parameters of the z ^ 7 LF given the characteristics of on- 
going surveys and models for potential future survey designs. 
The purview of our calculation was purposefully narrow since 
a more comprehensive evaluation of galaxy surveys could in- 
volve many additional questions we have not addressed. We 
now turn to a variety of possible caveats that stem from con- 
sidering photometric galaxy survey designs more generally. 

We have focused on forecasting constraints for the lu- 
minosity function. Our approach was modeled after 
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Fisher matrix calculations that used the abundance of 
galaxy c lusters to forec a st cosmological pa r amete r s con- 
straints (iHu & Kravtsov' '2003'; 'Lima^Hul l200l l2005t 
iCunha &Ev rard 2009; Wu et al. 2009), but other previous 
calculations have forecasted cosmological parameter con- 
straints from ga laxy clustering (e.g., Voaelev & Sz alav 1996,; 
Matsubara & Szalav 2001, 2003; Under 2003; Al brecht etal.l 
2009h . The incorporation of galaxy clustering data can cir- 
cumvent some assumptions made in ^when using simple 
abundance matching to assign galaxy bias by replacing the 
sample variance estimates in Equations [T5] andfTTjbv an inte- 
gral over the galaxy correlation function. Other estimates of 
how cosmic variance uncertainty is influenced by galaxy bias 
have taken a sknilar approach (e.g ., Newman & Davis 2002; 
Somerville et all 120041; IStark et al.ll2007bc iTrenti & Stiavellii 
20081k 

Our calculations have regarded a limited but interesting 
redshift regime near z ^ 7. While we have found that 
the combination of existing deep/narrow surveys with a fu- 
ture wide/shallow survey or a future ultradeep/narrow sur- 
vey would provide tight constraints on the z ^ 1 lumi- 
nosity function, studies of the galaxy population at higher 
and lower redshifts could require substantially different sur- 
veys. For instance, the z ^ 8 dropout candidates iden- 
tified in the UDF GO dat a are all fainter than //ab = 
27.7 (lBouwensetal.1 2009; B unker et al.ll2009l; iMcLureetal] 
l2OT0l:lYanetal.l 12009) . The decreasing abundance of rela- 
tively bright galaxies with increasing redshift will tend to fa- 
vor deeper and narrow surveys. Our calculations can easily be 
extended to estimate the constraining power of various sur- 
veys designs for higher-redshift galaxy populations, but we 
will save such estimates for future work when better fiducial 
estimates of the z > 8 luminosity function are available. 

Our Fisher matrix approach requires the use of a fiducial 
model for the abundance of z ^ 7. We adopt the Oesch et al] 
(1201 Obi) estimate of the galaxy luminosity function, which 
was determined by scaling the characteristic magnitude 
and normalization from lower redshift data and then fitting 
for the faint-end slope a. If the z ^ 7 galaxy luminosity func- 
tion differs substantially from the Oesch et al. (2010b) esti- 
mate, then our forecasted constraints could be similarly inac- 
curate. For instance, if the normalization (\>-^ was considerably 
lower or the ch aracteristic m agnitude much fainter than 
that estimate by lOesch et"aLl (l2010b ). then the relative bene- 
fit of combining flie UDF GO and GOODS ERS data with a 
wide/shallow survey over a narrow/ultradeep design could be 
reduced. 

The calculations in ^2.4.1| and ^suggest that splitting wide 
surveys into multiple fields to probe statistically-independent 
regions of the universe may not dramatically improve con- 
straints on the galaxy luminosity function. While this conclu- 
sion depends strongly on the total volume of the survey, other 
considerations such as scheduling, field observability, or sky 
backgrounds could make multiple fields advantageous com- 
pared with a single contiguous field of the same total area. 

The abundance matching calculation also requires either 
knowledge or assumption about the completeness of the sur- 
vey and the fraction of dark matter halos occupied by galax- 
ies. We have assumed that the surveys are essentially volume- 
limited and that there is a one-to-one correspondence between 
galaxies and dark matter halos. Both of these assumptions are 
likely imperfect, and some estimates of the high-redsh ift oc- 
cupation fraction are as low as 20% (IStark et al.ll2007bl) . The 
influence of these assumptions over the forecasted parame- 



ter constraints depends on the character of the survey. For 
a given observed luminosity function, reducing the halo oc- 
cupation fraction or the survey completeness acts to reduce 
the effective galaxy bias in the survey by either increasing 
the number of halos per galaxy in the survey or increasing 
the number of undetected galaxies. In either case, the Pois- 
son uncertainty is based on the observed number of galaxies 
and is unaffected. For purposes of constraining the observed 
luminosity function, wide surveys are fairly insensitive to ei- 
ther assumption, since Poisson uncertainty in the abundance 
of bright galaxies plays a large role in their error budget. The 
change in sample variance for narrow surveys can lead to a 
degradation of the marginalized parameter constraints (while 
the unmarginalized constraints can improve) by increasing 
the parameter correlations. However, the relative effect is 
small and the degradation is only 2 x for simultaneously low 
incompleteness (/comp = 0.1) and small occupation fraction 

(/occ=0.1). 

9. SUMMARY 

Motivated by the exciting initial galaxy survey data ob- 
tained by newly-installed Wide Field Camera 3 (WFC3) on 
the Hubble Space Telescope (HST), we have attempted to 
quantify how well on-going and possible future infrared sur- 
veys with WFC3 will constrain the abundance of galaxies at 
z ^ 7. Our primary methods and results include: 

• We perform a Fisher matrix calculation to forecast 
constraints on the galaxy luminosity function (LF) 
achievable by a survey with a given depth, area, and 
number of fields. In our approach, the constraints 
on the LF normalization 0*, characteristic magni- 
tude M^, and faint-end slope a that a survey can 
achieve directly relate to the sample cosmic variance 
and Poisson uncertainty on the observed galaxy abun- 
dance through the Fisher matrix. For a fiducial LF 
model, the abundance of observe d galaxies and dark 
matter halos are rnatche d (e.g., iConrov et al.l 120061; 
iConrov & Wechsleil2009l) to estimate the bias of galax- 
ies of a given luminosity. The galaxy bias is com- 
bined with the RMS density fluctuations within the sur- 
vey volume to calculate the sample cosmic variance 
(e.g., Newman & Davis 2002; Somerville et al. 2004 
IStark et al.ii2007bt iTrenti & Stiavellii i2008i) . while the 
Poisson variance simply scales with the square root of 
the number of of observed galaxies. The constraining 
power of each su rvey is then calcu lated from the Fisher 
matrix using the ISchechte^ ( 119761) model of the LF, its 
derivatives, the sample cosmic variance and Poisson 
uncertainties, and any data covariance. The combined 
constraints from multiple surveys can be estimated eas- 
ily by summing their Fisher matrices. Similar calcu- 
lations should prove useful for designing future photo- 
metric surveys and estimating their constraining power 
for the galaxy LF. 

• Using the Fisher matrix calculations, we estimate the 
constraints on the abundance of z ^ 7 galaxies that 
will be achieved with the entire forthcoming Ultradeep 
Field Guest Observation (UDF GO) and Great Obser- 
vatories Origins Deep Survey Early-Release Science 
(GOODS ERS) HST WFC3 IR channel data . Using 
the z ^ 7 galaxy LF estimated by Oesch et"an (l2010bb 
as a fiducial model, we calculate that the combined 
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UDF GO and GOODS ERS data will achieve marginal- 
ized (unmarginalized) LF parameter constraints of 
AM* « 0.5 mag (0.1 mag), Alogm^* « 0.4 (0.1), 
and Aa ~ 0.4 (0.15) when the surveys are fully com- 
pleted. These marginalized constraints correspond to 
uncertainties in the total number of z ^ 7 galaxies with 
magnitudes Muv < -18 (Muv < « -19.8) of 25% 
(200%), after accounting for covariances between the 
LF parameters. These surveys will provide the first 
detailed information on z ^ 7 galaxy populations, but 
abundance of the bright-end of the z ~ 7 LF will remain 
uncertain without further data. 



• We also forecast z ^ 7 LF constraints provided by a va- 
riety of model WFC3 surveys that would each require 
^ 450-900 HST orbits. The six model surveys con- 
sidered cover a large range of areas (14-3600 square 
arcmin) and depths (Hab =27-30 in F140W) to study 
the relative value area and depth for constraining the 
abundance of z ^ 7 galaxies. When combined with 
the forthcoming UDF GO and GOODS ERS data, all 
the surveys considered produce interesting luminosity 
function constraints (see Table |5]l. We find that a ^ 
1 deg^ survey to //ab ~ 27 in F140W provides the tight- 
est combined marginalized constraints (AM* w 0.14, 
A0* w 0.11, Aa « 0.20) on the abundance of z ^ 7 
galaxies of all survey designs we consider, but only by 
a small margin. This survey would require 780 total or- 
bits, including equal coverage in a bluer WFC3 filter to 
define a drop out color selection. In contrast, the abun- 
dance of faint galaxies would be best constrained by 
increasing depth of the HST ultradeep fields to ^ 125 
orbits per pointing (//ab ~ 30.1 in F140W), which 
provides marginalized LF constraints of AM* ss 0.22, 
A0* ss 0.15, and Aa w 0.10 for 750 total orbits (in- 
cluding equal coverage in a bluer WFC3 filter). 



• We also consider the usefulness of splitting surveys 
into A^eids multiple fields to probe independent sam- 
ples and re duce cosmic variance uncertainties (e.g., 
lNewman"& Davis 2002). We show that the shape of 
the ACDM power spectrum limits the statistical gain of 
splitting a high-redshift survey into multiple fields to 
< 10% (for A^fieids = 2) when the survey area is large (> 
0.5 deg^). We suggest that this statistical gain should be 
weighed against any scientific gains achieved by prob- 
ing large contiguous areas. 

Initial analyses of the UDF GO data have already demon- 
strated that the installation of WFC3 on HST will trans- 
form our knowledge of hi gh-redshift galaxy populatio ns 
at z ~ 7 and beyond (e.g., Bouwens et al. 2009, 2010a,'b; 
lOesch et al. 2010b a; Bunker et al. 2009; McLure et al. 2010; 
^iiet al. 2009; Wil kins et al.ll2010(: ILabbe et al.ll2009l 12010 : 
iFinkelste in et al. Our work has attempted to quantify 

expectations for the constraining power of the UDF GO and 
GOODS ERS surveys, and forecast constraints achievable 
with more extensive future surveys using WFC3 or other in- 
struments. These calculations illustrate how truly powerful 
the refurbished HST is for exploring high-redshift galaxy pop- 
ulations, and emphasize how exciting near-term gains in our 
knowledge of z > 7 galaxies will be. 
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Fig. 4. — Forecasted constr aints on z ' ■'.^ 7 lu minosity iiinction parameters expected from the existing GOODS ERS survey and the model Survey A. The 
constraints are calculated for a'Schechtei* 11910) luminosity function with a characteristic luminosity M*, normalization and faint-end slope a. Shown are the 
1 -cr constraints in the M» -logjQ (pi, (left panel), M^ — a (middle panel), and logjQ (f)-i, — a (right panel) space projections for the GOODS ERS (light blue region) 
and Survey A (blue region) surveys. Also shown are the constraints expected by combining both surveys (dark blue region). 
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Fig . 5 . — Forecasted constraints on z ~ 7 luminosity function parameters expected from the existing UDF GO survey and the model Survey B 1 . The constraints 
are calculated for a Schechter 1 1976) luminosity function with a characteristic luminosity M*, normalization c/)*, and faint-end slope a. Shown are the 1 - cr 
constraints in the M* - logjQ 0* (left panel), M* — a (middle panel), and log|Q(/>*-a (right panel) space projections for the UDF GO (light blue region) and 
Survey Bl (blue region) surveys. Also shown are the constraints expected by combining both surveys (dark blue region). 
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Fig. 6. — Forecast ed constraints on z ~ 7 luminosity function parameters expected from the existing UDF GO survey and the model Survey B2. The constraints 
ai'e calculated for a lSchechte3 119761) luminosity function with a characteristic luminosity M*, normalization and faint-end slope a. Shown are the 1 - cr 
constraints in the - logm 0* (left panel), M^ — a (middle panel), and logjQ(/>*-a (right panel) space projections for the UDF GO (light blue region) and 
Survey B2 (blue region) surveys. Also shown are the constraints expected by combining both surveys (dark blue region). 
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Fig. 7. — Forecast ed constraints on z ~ V luminosity function parameters expected from the existing UDF GO survey and tlie model Survey B3. The constraints 
are calculated for a lSchechte3 I197fl) luminosity function with a characteristic luminosity M*, normalization fA*, and faint-end slope a. Shown are the 1 -cr 
constraints in the M*-log]Q0* (left panel), M* — q (middle panel), and logm 0* - q (right panel) space projections for the UDF GO (light blue region) and 
Survey B3 (blue region) surveys. Also shown are the constraints expected by combining both surveys (dark blue region). 
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Fig. 8. — Forecasted constraints on z ~ 7 luminosity function parameters expected from the existing UDF GO and GOODS ERS surveys and the model 
Survey C. The constraints are calculated for a Schechter 1 1976) luminosity function with a characteristic luminosity M,, normalization (f)^, and faint-end slope 
a. Shown are the l-cr constraints in the M* — logjQ^* (left panel), M„-a (middle panel), and logjQflii — a (right panel) space projections for the combined 
UDF GO and GOODS ERS surveys (light blue region) and the model Survey C (blue region). Also shown are the constraints expected by combining all three 
surveys (dai'k blue region). 
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Fig. 9. — Forecasted constraints on z ~ 7 luminosity function parameters expected from the existing UDF GO and GOODS ERS surveys and the model 
Survey D. The constraints are calculated for a Schechter ( 1976) luminosity function with a characteristic luminosity M*, normalization 0*, and faint-end slope 
Q. Shown are the I - a constraints in the M* - logjQ 0* (left panel), - a (middle panel), and log[Q 4'*~oi (right panel) space projections for the combined 
UDF GO and GOODS ERS surveys (light blue region) and the model Survey D (blue region). Also shown are the constraints expected by combining all three 
surveys (dark blue region). 
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TABLE 4 

Forecasted Constraints for Existing Surveys 



Survey 


Marg. Uncert. 


Unmarg. Uncert. 


Marg. Uncert. 


Unmarg. Uncert. 


Marg. Uncert. 


Unmarg. Uncert. 


Name 


Mt, 


M« 


logio "P* 


logio <P* 


a 


a 


UDF GO 


0.837 


0.160 


0.655 


0.093 


0.566 


0.162 


GOODS ERS 


1.338 


0.171 


1.026 


0.145 


1.852 


0.520 


Combined 


0.524 


0.117 


0.425 


0.078 


0.415 


0.154 


Survey 


Pearson p 


Pearson p 


Pearson p 


&(M^,-(j>^,) 


0(M^,-a) 


©((^* - a) 


Name 


M*-logio</>* 


Mi,-a 


logio 9^*-" 


[deg] 


[deg] 


[deg] 


UDF GO 


0.98 


0.91 


0.95 


38 


33 


41 


GOODS ERS 


0.99 


0.96 


0.95 


37 


55 


62 


Combined 


0.97 


0.89 


0.93 


39 


38 


44 








TABLE 5 












Forecasted Constraints for Model Surveys 






— 

Survey 


Marg. Uncert. 


Unmarg. Uncert. 


Marg. Uncert. 


Unmarg. Uncert. 


Marg. Uncert. 


Unmarg. Uncert. 


Name 


)i/f 


M^, 


logio </>* 


logio 0* 


CK 


Q; 


Survey A 


0.276 


0.158 


0.163 


0.064 


0.104 


0.057 


Survey B 1 


0.225 


0.084 


0.185 


0.053 


0.182 


0.055 


Survey B2 


0.257 


0.084 


0.210 


0.057 


0.231 


0.068 


Survey B3 


0.299 


0.085 


0.244 


0.061 


0.305 


0.089 


Survey C 


0.244 


0.052 


0.193 


0.041 


0.311 


0.088 


Survey D 


0.336 


0.041 


0.194 


0.045 


0.722 


0.137 


Ex.+Survey A" 


0.221 


0.116 


0.151 


0.058 


0.101 


0.057 


Ex.+Survey Bl* 


0.214 


0.074 


0.176 


0.046 


0.172 


0.052 


Ex.+Survey Bl* 


0.236 


0.074 


0.194 


0.049 


0.207 


0.063 


Ex.+Survey BS* 


0.253 


0.075 


0.210 


0.051 


0.245 


0.078 


Ex.+Survey 


0.185 


0.048 


0.150 


0.036 


0.219 


0.077 


Ex.+Survey D*^ 


0.136 


0.039 


0.111 


0.039 


0.201 


0.103 


Siirvpv 




Pearson p 


Pearson p 








Name 


M* -log in (b* 




logio </"*"" 


[deg] 


[deg] 


[deg] 


Survey A 


0.79 


0.52 


0.81 


28 


12 


30 


Survey B 1 


0.92 


0.91 


0.95 


39 


38 


45 


Survey B2 


0.94 


0.93 


0.95 


39 


42 


48 


Survey B3 


0.95 


0.94 


0.95 


39 


46 


52 


Survey C 


0.97 


0.95 


0.95 


38 


52 


59 


Survey D 


0.96 


0.97 


0.90 


29 


65 


76 


Ex.+Survey A" 


0.84 


0.62 


0.82 


33 


18 


31 


Ex.+Survey Bl'' 


0.94 


0.91 


0.95 


39 


38 


44 


Ex.+Survey B2'' 


0.95 


0.92 


0.95 


39 


41 


47 


Ex.+Survey B3'' 


0.95 


0.92 


0.95 


39 


44 


50 


Ex.+Survey 


0.96 


0.92 


0.93 


39 


50 


56 


Ex.+Survey D'^ 


0.93 


0.86 


0.77 


39 


58 


65 



" Combined with the existing GOODS ERS survey 

Combined with the existing UDF GO survey. 
<^ Combined with the existing UDF GO and GOODS ERS surveys. 



